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ABSTRACT 

Velocity gradients in a stellar atmospheric plasma have an effect on the anisotropy of the radiation 
field that illuminates each point within the medium, and this may in principle influence the scattering 
line polarization that results from the induced atomic level polarization. Here we analyze the emergent 
linear polarization profiles of the Ca n infrared triplet after solving the radiative transfer problem of 
scattering polarization in time-dependent hydrodynamical models of the solar chromosphere, taking 
into account the effect of the plasma macroscopic velocity on the atomic level polarization. We discuss 
the influence that the velocity and temperature shocks in the considered chromospheric models have 
on the temporal evolution of the scattering polarization signals of the Ca n infrared lines, as well 
as on the temporally averaged profiles. Our results indicate that the increase of the linear polariza- 
tion amplitudes caused by macroscopic velocity gradients may be significant in realistic situations. 
We also study the effect of the integration time, the micro-turbulent velocity and the photosphcric 
dynamical conditions, and discuss the feasibility of observing with large-aperture telescopes the tem- 
poral variation of the scattering polarization profiles. Finally, we explore the possibility of using the 
differential Hanle effect in the IR triplet of Ca n with the intention of avoiding the characterization 
of the zero-field polarization to infer magnetic fields in dynamic situations. 

Subject headings: Polarization - scattering - radiative transfer Sun: chromosphere Stars: atmospheres, 
shocks 



1. INTRODUCTION 

The chromosphere, the interface region between the 
photosphere and the corona, is a very important part 
of the solar atmosphere. It is the place where most 
of the non-thermal energy that creates the corona and 
solar wind is released, with a heating rate requirement 
that is between one and two orders of magnitude larger 
than in the corona. To infer the thermal, dynamic and 
magnetic structure of the solar chromosphere is thus a 
very important goal in astrophysics. For instance, it 
is believed that the dissipation of magnetic energy in 
the 10 6 K corona may be significantly modulated by the 
strength and st ructure of th e magnetic field in the chro- 
mosphere (e.g., Parker 2007). However, "measuring" the 
chromosph eric magnetic field is notoriously difficult (e.g., 
reviews by ICasini fc Landi Degrinno ccnti 2002 iHarvevi 
I2009t ITruiillo Buenoll2010D . While spectroscopic obser- 
vations allow us to determine temperatures, flows and 
waves, they do not provide any quantitative informa- 
tion on the chromospheric magnetic field. To this end, 
we need to measure and interpret the polarization t hat 
some physical mechanisms introduce in chromospheric 
spectral lines. These mechanisms are the Zeeman effect, 
scattering processes and the Hanle effect. 

The circular and linear polarization signals that the 
Zeeman effect can in principle produce in a spectral 
line are caused by the wavelength shifts between the 
7r and a transitions of the line, as a result of the Zee- 
man splitting induced by the presence of a magnetic 
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field. The amplitude of the circular polarization scales 
with the ratio, 1Z, between the Zeeman splitting and the 
Doppler line width. Th e amplitude of the linear polariza- 
tion s cales with TZ 2 (see lLandi DeglTnnocenti fc Landolfil 
l2004f ). Outside sunspots (where i?<100 G at chromo- 
spheric heights) which explains why it is so dif- 
ficult to detect the polarization of the Zeeman effect in 
a chromospheric line. Typically, only the circular polar- 
ization is detected, especially in long- wavelength chromo- 
spheric lines such as those of the IR triplet of Ca II (e.g., 
ITruiillo Bueno|[2010L figure 3). But the linear polariza- 
tion observed in quiet regions of the solar chromosphere 
has practically nothing to do with the transverse Zeeman 
effect. 

In weakly magnetized regions, the linear polarization 
of chromospheric lines is dominated by scattering pro- 
cesses. The physical origin of this polarization is the 
difference among the electronic populations of sublevels 
pertaining to the levels of the spectral line under consid- 
eration. This so-called atomic level polarization, which is 
caused by the anisotropic illumination of the atoms, pro- 
duces selective emission and/or selective absorption of 
polarizatio n components without the need of a magneti c 
field (e.g.. iManso Sainz fc Truiillo Buenol l2003bl . 120101) . 
The larger the anisotropy of the incident radiation field 
the larger the induced atomic level polarization and the 
larger the amplitude of the linear polarization of the 
emergent spectral line radiation. In an optically thick 
plasma like the solar atmosphere, the anisotropy of the 
radiation field depends mainly on the spatial distribu- 
tion of the physical quantities that determine, at each 
point within the medium, the angular variation of the 
incident intensity. Great attention has been paid to the 
gradient of the source function (e g..lTruiillo Buendl200TI : 
lLandi Pegl'Innocenti fc Landolfil 120041 ) but, in a highly 
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dynamic medium like the solar chromosphere, the gra- 
dients of the macroscopic v elocity of the plas ma may 
also play an important rol e dCarlin et al1l2012l and ref- 
erences therein). In fact, in lCarlin et al.l (|2012L hereafter 
Paper i) we showed that it can affect significantly the 
scattering polarization of the IR triplet of Ca II. Our 
arguments were based on radiative transfer calculations 
in a semi-empirical model of the solar atmosphere, af- 
ter introducing ad-hoc velocity gradients and comparing 
the computed Q/I profiles with those corresponding to 
the static case. Given the diagnostic potential of the 
Ca II IR triplet for exploring the magnetism of the solar 
chromosphere (e.g..lManso Sainz fc Truiillo Buendl20ll 
iDe la Cruz Rodriguez et al.ll2012f ) . and the fact that the 
region where such chromosphcric lines originate may be 
affected by vigorous a nd repetitive shock waves (e.g., 
iCarlsson fc Steir3ll997[ ) , it is necessary to investigate the 
radiative trans fcr problem of scattering polarization in 
the Ca II IR triplet using dynamical, time-dependent at- 
mospheric models of the solar chromosphere. In this pa- 
per, we show the results of such investigation. 

2. DESCRIPTION OF THE PROBLEM AND THE 
RESOLUTION PROCEDURE. 

We have carried out radiative transfer calculations of 
the linear polarization produced by scattering in the 
Ca II infrared (IR) triplet. The polarization is pro- 
duced by the atomic level polarization that results from 
anisotropic radiation pumping in the hydrodynamical 
( HP) models of sola r chromosph eric dynamics described 
in ICarlsson fc Steinl ffl997L l2T)f)l . 

We used two time series of snapshots from the above- 
mentioned radiation HD simulations, each one lasting 
about 3600 s and showing the upward propagation of 
acoustic wave trains growing in amplitude with height 
until they eventually produce shocks. The first one cor- 
responds to a relatively strong photosphcric disturbance 
showing well-developed cool phases a nd pronounced hot 
zones at chromosphcric heights (see ICarlsson fc Steinl 
Il997t we refer to this as the strongly dynamic case). The 
second simulation corre sponds to a less intense photo- 
spheric disturbance (see ICarlsson fc Steinl . I2002t we re- 
fer to this as the weakly dynamic case). Thus, the ther- 
modynamical evolution of the atmosphere (including the 
chromosphere and the transition region) is driven by the 
bottom boundary condition that is imposed on the veloc- 
ity. This realistic boundary condition is extracted from 
the measured Doppler shifts in the Fc I line at 3966.8 A. 
Our description focuses mainly on the strongly dynamic 
case, but in Sec. 14.51 wc compare the results with those 
corresponding to the weakly dynamic case. 

To characterize the simulations we can use the follow- 
ing quantities. In terms of the velocity gradients, and 
using units related to a representative scale height □ V. = 
275 km, the temporal average of the maximum velocity 
gradient along the atmosphere is 40 km • s _1 per scale 
height (or 145 m • s _1 km _1 ) in the strongly dynamic case 

4 A scale height can be denned as the typical distance over which 
atmospheric magnitudes such as the density vary an order ol mag- 
nitude. But the scale height is not a fixed quantity in a time- 
dependent model that contains important temporal variations in 
those magnitudes. For this reason we have defined an averaged 
scale height as the representative value used for the characteriza- 
tion of the velocity gradients. 



and 13km-s _1 per scale height (or 47m-s~ 1 km~ 1 ) in 
the weakly dynamic case. Likewise, the temporal aver- 
age of the minimum of temperature in the atmosphere is 
3976 K in the strongly dynamic case, and 4292 K in the 
weakly dynamic case. 

At each time step of the HD simulation under consid- 
eration we use the corresponding one-dimensional strati- 
fications of the vertical velocity, temperature and density 
to compute the emergent J(A) and Q(X) profiles through 
the application of the multilevel r a diative transf er code 
of IManso Sainz fc Truiillo Buenol (|2003al l20Tot) . after 
the generali z ation to the non-static case described in 
iCarlin et all (|2012ft . Specifically, we have solved jointly 
the radiative transfer (RT) equations for the Stokes / and 
Q parameters and the statistical equilibrium equations 
(SEE) for the atomic populations of each energy level 
and the population imbalances among its magnetic en- 
ergy sublevels (equivalently, the multipolar tensor com- 
ponents of the atomic density matrix, p^(Ji), with Jj 
the angular momentum of each level i). This is the 
NLTE radiative transfe r problem of the second kind (see 
sectio ns 7.2 and 7.13 in lLandi Degl'Innocenti fc Landoffil 
200 1 : :. Once the self-consistent solution of such equa- 
tions is found at each height in the atmospheric model 
under consideration, we compute the coefficients of the 
emission vector and of the propagation matrix (see sec- 
tion 2.2 of Paper I) and solve the RT equations for a 
line of sight (LOS) with /i = 0.1, where fi is the cosine 
of the heliocentric angle. This LOS has been chosen in 
order to simulate a close to the limb obse r vation , such as 
that shown in figure 13 of iStenflo et al.l (|2000[) . To ac- 
count for the macroscopic motions, we have introduced 
the Doppler effect in the calculation of the absorption 
and emission profiles for each wavelength and ray direc- 
tion (Paper i). The influence of the Doppler effect on the 
SEE appears directly because the radiative rates depend 
on the radiation field tensor components. Likewise, the 
RTE is affected because the Doppler effect modifies the 
elements of the propagation matrix and of the emission 
vector. 

Given that the computations reported here are carried 
out in plane-parallel atmospheric models, it is necessary 
to introduce a micro-turbulent velocity, that accounts for 
the Doppler shifts (inducing an effective line broadening) 
produced by moving fluid elements below the resolution 
element. In order to estimate a suitable value (assumed 
constant with height), we have calculated the emergent 
intensities at disk center and compared th em with those 
of the solar Kitt Peak FTS Spectral Atlas (|Kurucz et al.l 
I1984D . A good agreement is obtained with 3.5 km s _1 . 

3. DESCRIPTION AND CHARACTERIZATION OF THE 
RESULTS. 

A standard Fourier analysis of the atmosphere model 
shows that it acts as a pass-band filter for the multifrc- 
quency sound waves generated in the lower boundary. 
The result is that the predominating periods at chro- 
mospheric heights and h igher are around three minutes 
(jCarlsson fc Steinl [1997) . For practical reasons we di- 
vided the temporal evolution in 3-min intervals so that 
the beginning of each interval coincides with the moment 
in which the shock front in temperature and velocity is 
the sharpest in each interval (vertical lines in figures with 
temporal axis, like Fig. [1]). Given the power of the 3-min 
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waves, this division turns out to be "natural" and can be 
used to mark the most interesting events we see in the 
emergent polarization. 

Inside each three-minutes cycle we distinguish between 
compression and expansion phases. They can be easily 
identified following the height at which t 1 ° q s = 1, i.e., 
where the optical depth at line center (v ) along the LOS 
equals unity (upper panel of Fig. []}. This quantity is a 
good marker of the shock fronts when they cross heights 
between 1 and 2 Mm. It is because the steep changes 
in opacity inside the shocks forces the r = 1 region to 
remain comprised within them. The line transitions at 
8542 A and 8662 A (green and red lines in the upper 
panel of Fig. [I} follow a clearer periodic pattern because 
they form higher, where less frequency components of 
the velocity waves can arrive. Compression phases be- 
gin when plasma falls down from upper layers (rg^ = 1 
and Tggg 2 = 1 decrease in top panel of Fig. []}, while 
simultaneously a new upward propagating wave emerges 
amplified into the chromosphere. At the end of this stage 
a shock wave is completely developed and the rl° s = 1 
position is close to ~1200 km for the three IR lines. The 
shock waves so created start always in such region be- 
tween 1 and 1.5 Mm0. After that, during what we term 
expansion phase (heights for Tgg| 2 = 1 and Tggg 2 = 1 aris- 
ing in top panel of Fig. [IJ, the shock fronts travel upward 
increasing the plasma velocities as they encounter lower 
densities. 

Figure Q] also shows the time evolution of other quan- 
tities during the first 2000 s after the initial transient. 
In the second row, the location and value of the tem- 
perature minimum are displayed, showing a clear corre- 
spondence with expansion and contraction phases. In 
the third row, we show the ensuing variation of (Q/I) pp , 
defined as the peak-to-peak difference of the Q/I pro- 
file for each spectral line. It is a measure of the lin- 
ear polarization signal contrast that was used in Paper I 
to characterize the polarization amplitude and discrim- 
inate their variations with respect to static cases. In 
each cycle we see an amplification of (Q/I) pp occurring 
at expansion phases and an usually larger amplification 
during contraction phases. Finally, the time evolution 
of the emergent Q(A)/I(A) profile for the 8542 A line 
is illustrated in the lower panel (the vertical axis shows 
0.6 A around the rest wavelength of the line). Here, 
we observe two distinct areas showing amplifications in- 
side each three-minute cycle. The first amplification is 
blue-shifted, because it happens in an atmospheric ex- 
pansion phase (plasma moving towards the observer) . It 
is weaker than the second amplification, which is red- 
shifted and occurs during the compression phase (plasma 
moving down in the atmosphere) . This indicates that the 
compression phase is more efficient producing a polariza- 
tion amplification than the expansion one. The reason is 
that during compression we have stronger velocity and 
temperature gradients along the main regions of forma- 
tion. Following the results of Paper I, the larger the 
gradient, the larger the enhancement of the linear po- 
larization signal. The behaviour is similar in the other 
transitions. 

5 It is in this range of heights where the Ca II IR triplet forms 
in typical semi-empirical models 



There is a clear correspondence between the maximum 
value of the temperature minimum (hot-chromosphcre 
time-steps) and the largest peaks of the (Q/I) pp signal, 
taking place just before the maximum contraction (dot- 
ted vertical lines). As the atmosphere is compressed, 
the temperature increases at chromospheric heights and 
the resulting gradient of the source function produces 
an increase of the radiation field anisotropy in the up- 
per layers. This directly leads to an enhanced emer- 
gent linear polarization signal. On the contrary, in cold- 
chromosphere models the expansion reaches its maxi- 
mum and (Q/I) PP is near its minimum value. 

Even in such complex situations, we still witness the 
already known effects of amplification (with respect to 
the static case), frequency shift and asymmetry in the 
linear polarization profiles due to dynamics. All of them 
have been already explained in Paper I, using the semi- 
empirical FAL-C model of Fontenla et al. (1993) with 
ad-hoc velocity stratifications. The enhancement is pro- 
duced as a consequence of the velocity gradients and sub- 
sequent anisotropy enhancements. However, some differ- 
ences exist from the experiments in semi-empirical mod- 
els and the calculations presented in this paper. First, 
the velocity stratification in the HD models is, in gen- 
eral, non-monotonic and with a non-constant variation 
with height. Second, the maximum velocity gradients 
are located at shocks, with amplitudes that reach tens or 
even hundreds of meters per second per kilometer (as a 
comparison, in Paper I we dealt with velocity gradients 
between and 20 m s _1 km -1 ). Third, as commented be- 
fore, we have shocks in temperature that produce larger 
source function gradients and additional enhancement of 
the radiation anisotropy and of the linear polarization. 
Finally, these variations are usually concentrated in the 
formation regions of the triplet lines. All these mecha- 
nisms act together and enhance the linear polarization 
of the emergent radiation with amplification factors up 
to - 10 (in the 8498 Aline) and ~ 7 (in the 8542 A and 
8662 A lines), for the instantaneous values of the Q/I 
amplitudes with respect to the static FAL-C case. How- 
ever, if we consider temporal averages of the emergent 
Stokes profiles during long periods, we get amplification 
factors of about a factor of 2 (time- averaged Q/I ampli- 
tudes reach - 1 % for 8542 A and 8662 A lines). 

Summarizing, the temporal evolution of the polariza- 
tion is driven by the temperature and velocity stratifica- 
tions, that in turn are a result of the dynamical condi- 
tions set in the photosphere. 

4. ANALYSIS AND DISCUSSION OF RESULTS. 

4.1. The effect of the velocity. 

A way of visualizing the effect of vertical velocity gra- 
dients on the emergent scattering polarization is to com- 
pare the evolution of the polarization profiles correspond- 
ing to both the static and non-static case. In the ab- 
sence of velocities (lower row of Fig. [2]), the maximum 
of the Q/I profiles is always located at A = A (i.e., line 
center), and its temporal evolution presents a sawtooth 
shape. When the effect of velocities is included in the 
calculations (upper row of Fig. [5]), the maximum of the 
Q/I signal is no longer located at the central wavelength 
and its temporal evolution assumes a different shape with 
two peaks every 3-min period (upper right panel). These 
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Fig. 1. — Top row: time evolution of the atmospheric heights where T l ° a = 1 for the three Ca II IR transitions. The dotted vertical 

lines are located at the local minimum of the Tg^j = 1 curve, and they can be considered to indicate the beginning and the end of each 
"three-minute" period with their corresponding expansion and compression phases. Second row: time evolution of the temperature value 
and atmospheric height of the temperature minimum. Third row: time evolution of the polarization contrast (max(Q/I) — min(Q/I)) for 
the three lines of the Ca II IR triplet. The polarization amplitude of the A8498 line has been multiplied by -5 to show the results for the 
three lines on the same scale. Note that, by definition, the line contrast is always positive, so with this quantity we cannot know if the 
polarization signal is positive or negative. Thus, with the artificial sign inversion done for the 8498 A line contrast in this figure, we try 
to illustrate that this line is the only one whose larger polarization amplitudes are negative. Bottom row: time evolution of the calculated 
Q(X)/ I(X) fractional linear polarization profile of the A8542 line. 

of the velocity gradients in the HD models, we see that 
the resulting temporally-averaged scattering polarization 
signals (which include the impact of the temperature and 
density shocks) arc similar to the Q/I profiles computed 
in the static FAL-C semi-empirical model. 

4.2. The combined effect of velocity and temperature on 
the linear polarization. 

In Fig. [3] we display some relevant magnitudes for 
three different situations in the simulation. The first col- 
umn corresponds to a quiet time-step, with no shocks, 
zero velocity and without any kind of amplifications (it 



wavelength and amplitude modulations are produced by 
the Doppler effect of the velocity gradients. 

It is interesting to compare the mean Q/I amplitudes 
obtained in the hydrodynamical models with the one cal- 
culated in the FAL-C model. They differ notably (see 
horizontal lines in right panels of Fig. [5]). In the 8542 A 
transition we have mean values around 1%, 0.31% and 
0.42% for the HD models with velocities, the HD mod- 
els at rest and static FALC models, respectively. The 
results of these figures have been obtained using an inte- 
gration time of 1040 s (^ 17 minutes), the duration of the 
temporal interval shown in Fig. [5J Neglecting the effect 
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Fig. 2. — Temporal evolution of Q(A)//(A) (left column panels) and of the Max(Q/7) peak amplitude (right column panels) during 
1040 s (17 minutes) for the 8542 A transition. The results of the upper panels take into account the effect of the velocity gradients, 
while those of the lower panels were obtained assuming zero macroscopic velocities while calculating the density-matrix elements. The 
vertical lines indicate the beginning and the end of each three-minutes period. The horizontal lines in the right column panels show the 
temporally-averaged amplitudes {(Q) / (I)) obtained from the Q/I profiles computed in the strongly dynamic model, with the solid line 
indicating the polarization amplitude obtained when the effect of the velocity gradients is taken into account and the dashed lines showing 
the polarization amplitude obtained assuming zero macroscopic velocities. The dotted horizontal lines in the right panels indicate the 
polarization amplitudes obtained in the (static) FAL-C semi-empirical atmospheric model. 

the source function. Additionally, during contraction the 
temperature shock occurs in optically thick and denser 
layers (deeper layers below t 1 ° s = 1), forcing the source 
function gradient to increase with respect to the static 
case at those heights. Note how in this last case the 
source function rises as a whole because of the warm- 
ing (compare the source function in the middle panel, 
the black solid line that has been obtained neglecting 
velocities, with the non-dynamic source function in the 
left column). If the macroscopic velocity is now consid- 
ered, we additionally get a jump in the source function 
(red lines in middle column of Fig. [3]) caused by the ve- 
locity shock that is developed in this contraction phase. 
This behaviour is accompanied by a significant Dopplcr- 
induccd anisotropy enhancement that amplifies the linear 
polarization, as shown in the corresponding lower panels 
of the same figure. 

In the expansion phases, the shock waves move upward 
and the chromosphere becomes cooler. This induces a 
lower source function and smaller polarization ampli- 
tudes (as compared with the contraction phase). Oth- 
erwise, as the density of scatterers is now lower around 
the shock (because it moved upward to regions with 
tI° s < 1), the temperature gradients have smaller effects 
on the polarization profiles than during the contraction 
phases. In this expansion time step, the black solid line 
representing the static source function is similar to the 
non-dynamic source function of the left column. How- 
ever, once the motions are introduced, and despite of the 
fact that the shocks have already reached upper chromo- 
sphcric layers, the remanent velocity field has still a siz- 



is the initial transient phase). The central column shows 
a phase of compression, in which shocks are important. 
Finally, the last column displays an expansion phase, in 
which the atmosphere is expanded and the shocks are al- 
ready travelling over the transition region. Furthermore, 
we distinguish between the solutions when motions are 
taken into account (red lines) and the solutions obtained 
allowing shocks in all magnitudes but artificially setting 
the velocity to zero (black lines). 

The normalized velocity £ z = (vq/c)v z / Auo-, with 
A^d the Dopplcr width of the absorption profiles (that 
depends on the temperature), c the speed of light and v z 
the vertical velocity, is the quantity that controls the im- 
portance of the atmospheric motions in relation with the 
radiation anisotropy and the scattering polarization (see 
Paper i). Note that this quantity considers the combined 
effect of velocity and temperature. In the HD atmosphere 
models, £ z (solid lines in upper panels of Fig. [3J is only 
significant in the formation region of the IR triplet lines 
( T l° s ~ 1 region with high velocity gradients, and not 
very high temperatures). Although shock waves increase 
the chromosphcric temperature, the effect of the velocity 
gradients is predominant. The contrary occurs over the 
transition region, where the thermal line width is much 
larger than the Doppler shifts. 

The expansion and contraction can be identified also in 
quantities such as the intensity source function and the 
Planck function (second row in Fig. [3]). During contrac- 
tion phases (middle column panels), high temperatures 
produce a more efficient population pumping towards up- 
per levels, incrementing the emissivity and, consequently, 
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Fig. 3. — This figure shows the values of the quantities indicated in the vertical axes for the time steps 10 s. (left column, a 'quiet' situation). 
770 s. (middle column, a compression stage) and 870 s. (right column, an expansion phase) of the strongly dynamic hydrodynamical 
simulation. First row (from top): macroscopic vertical velocity (dotted line) and adimcnsional vertical velocity £ z (solid line) versus height. 
The green line of each panel indicates the atmospheric height where t^° s = 1. Second row (from top): temperature (dashed line), Planck 
function (dotted line), source function for the zero velocity approximation (solid black line) and source function allowing the influence of 
the model's velocity gradients (red solid line). The green lines mark the instantaneous positions of tI° s = 1 and rl° s = 1. Third row (from 
top): line anisotropy factor (see equation 4 in Carlin et al. 2012) calculated for each of the above-mentioned cases; neglecting the effect 
of the velocity gradients (black solid line) or allowing it (red solid lines). Fourth row (from top): emergent Q/I profiles versus wavelength 
with respect to the line center, with the same color code than in previous panels. The dashed lines here are the time averages over the 
entire simulation, for each case (with or without velocities). 

able gradient that enhances the source function (Doppler 4.3. Averaged values of the polarization profiles. 

brightening effect). 
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Fig. 4. — Fractional linear polarization profiles of the 8498 A 
and 8542 A lines after temporally averaging the Stokes I and Q 
profiles during 3070 seconds (51 minutes). These (Q)/(J) pro- 
files may be considered to emulate what can be actually observed 
with today's solar telescopes. Black solid profiles: static case with 
'-'micro = 3.5kms~ 1 . Red solid profiles: strongly dynamic case 
taking into account the effect of the velocity gradients and as- 
suming ti m i cr o = 3.5kms _1 . Black dashed profiles: strongly dy- 
namic case neglecting the effect of the velocity gradients and as- 
suming ?) m icro = O.Red dashed profiles: strongly dynamic case 
taking into account the effect of the velocity gradients and assum- 
ing i> m i C ro = 0. The green solid lines show the temporally aver- 
aged profiles obtained after applying the velocity free approxima- 
tion (VFA) with 

^micro — 3.5 km s (i.e., neglecting the Dopplcr 
shifts of the macroscopic velocities when computing the density 
matrix elements, but taking them into account when calculating 
the emergent Stokes profiles). 

In order to compute the average linear polarization sig- 
nal that one would observe without any temporal resolu- 
tion, we average Q and I (obtaining (Q) / (I)) over 3070 s 
(« 51 minutes) for four different cases (Fig. We con- 
sider the cases with zero micro-turbulent velocity (dotted 
lines) and a constant micro-turbulent velocity of 3.5 km 
s _1 (solid lines). For each case, we distinguish between 
the results switching off the velocity (black lines) and 
the results allowing for macroscopic velocity fields (red 
lines) . 

When macroscopic motions are considered, the polar- 
ization profiles become asymmetric. Furthermore, they 
become more negative in the case of the 8498 A transi- 
tion and more positive in the other two transitions. The 
asymmetry of the red profiles is a consequence of the fact 
that, during the averaging period, the dynamical situa- 
tions in which the velocity gradient is negative (veloc- 
ity field mostly decreasing with height) dominate over 
the situations with velocity gradients that are mostly 



positive. This predominance is not because the situa- 
tions with negative velocity gradients are more frequent 
but because such situations are more efficient on ampli- 
fying the linear polarization. This happens during the 
compression phase because i) the velocity gradients are 
larger, ii) there is also a shock in temperature affecting 
the formation region and iii) the shock fronts are located 
just below the t 1 ° q s = 1 height. The results are quali- 
tatively the same independently of the micro-turbulent 
velocity value, but, when it is not considered, the ampli- 
fication of (Q/I) pp is larger and the profiles are narrower. 

If we decrease the averaging interval to 9 minutes, we 
obtain profiles that are essentially similar to the ones ob- 
tained by averaging during 51 minutes (showed in Fig. 

If we integrate less than that, significant variations 
appear in the shape and amplitude of the emergent pro- 
files. This indicates that, concerning the linear polariza- 
tion, there is still reliable dynamic information contained 
in a time interval corresponding to a few 3-min cycles. 

4.4. The velocity free approximation. 

An approximation that is sometimes applied to solve 
radiative transfer problems in dynamical atmospheres 
(either taking into account the presence of atomic polar- 
ization or not) is the velocity free approximation (VFA). 
It is based on solving the SEE and RTE simultaneously 
but neglecting the effect of plasma motions. However, 
once they are consistently solved, such plasma motions 
are included in the synthesis of the emergent Stokes pro- 
files (along [i — 0.1 in our case). Consequently, the 
density matrix elements are calculated as if plasma mo- 
tions did not affect them, reducing the complexity and 
computational effort of the problem since a reduced fre- 
quency grid is used to compute the mean intensity and 
the anisotropy. The results of applying it to each time 
step of our HD evolution is the temporal average illus- 
trated as the green line in Fig. 2] This approximation is 
clearly not appropriate in our case, given that the pro- 
files just become asymmetric (with respect to the static 
profiles) but without the amplification. The reason for 
this lack of amplification is that the anisotropy control- 
ling the linear polarization is not correctly enhanced (sec 
Paper i). On the other hand, the asymmetry is purely 
due to the asymmetric absorption with respect to the line 
center that motions produce along the ray under consid- 
eration. Hence, in order to obtain reliable results it is 
mandatory to include the effect of Doppler shifts in the 
whole set of equations, and we conclude that the VFA 
should not be applied. 

4.5. The effect of photospheric dynamics. 

Given that the small velocity fields appearing in the 
photosphere are amplified because of the exponential de- 
crease in the density while the perturbations travel out- 
wards, the properties of the bottom boundary condi- 
tion are determinant in the behaviour of the emergent 
Stokes parameter of chromosphcric lines. We compare 
the strongly dynamic case that forms the core of our pa- 
per with the weakly dynamic case that has been already 
introduced in Sec. [3] Although the mean maximum 
velocity gradient is three times smaller in the weakly 
dynamic case and the averaged polarization amplitudes 
are also smaller than in the strongly dynamic one, we 
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Fig. 6. — Comparison of the results in the strongly and weakly dynamic cases (see Sec. 14.51 1. All profiles are the resulting (Q) / (I) profiles 
obtained averaging Q and / during 15 minutes of the considered simulation. Solid lines: results calculated in the strongly dynamic case. 
Dashed lines: results calculated in the weakly dynamic case. Black lines: results allowing variations in all magnitudes but neglecting the 
velocity. Red lines: results allowing the effect of the velocity gradients. 
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Fig. 7. — Temporal evolution of the emergent 8542A Q/I profiles for different integration times. Wavelength is in the vertical axis. From 
top to bottom we have a 10, 60 and 180 s. temporal resolution, respectively. The white solid lines show the temporal evolution of (Q/I) PP 
(i.e., the amplitude contrast at each time-step). Vertical dotted lines mark each three-minute period. 



still find comparable or even slightly larger instantaneous 
(Q/I) pp amplitudes (see Fig. [5]). The resulting averaged 
polarization profiles arc qualitatively the same but they 
differ in amplitude (Fig. [6|). This is a reasonable re- 
sult because in the weakly dynamic scenario the instan- 
taneous velocity gradients are smaller in general. Dif- 
ferences are especially critical for the 8498 A line, whose 
linear polarization profiles can be positive, but also adopt 
significant negative values at redder wavelengths (when 
velocity gradients are mainly positive with height) or at 
bluer wavelengths (when velocity gradients are mainly 
negative with height). This behaviour produces cancel- 



lation effects with integration times larger than a three- 
minute period. Furthermore, the central depression pro- 
duced in the 8498 A average profile when the velocity 
is neglected in the strongly dynamic simulation (solid 
black line in left panel of Fig. O do not appear in the 
average profile corresponding to the weakly dynamic case 
(dashed black line in the same panel) because of the dif- 
ferences in the instantaneous temperature stratifications. 
The sensitivity of this spectral line to the instantaneous 
photospheric perturbations and to the developed chro- 
mospheric shocks is larger than in the other two lines. 
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Fig. 8. — Resulting Q/I profiles after averaging Q and I during 15 
minutes in the strongly dynamic simulation, using different values 
of F. F is the scaling factor by which we have multiplied the 
modulus of the macroscopic velocity at each atmospheric height. 
The curves correspond to F = 1, 0.9, 0.7, 0.5, 0.2, 0.1, 0, going 
from black colour (fully dynamic case) to lighter grey colour (static 
case). The results for the 8662 A line are very similar to the ones 
obtained for the 8542 A line. 

4.6. The effect of the integration time. 

In order to detect in the Sun the time evolution of the 
linear polarization signals, the observations must have 
enough time resolution, signal to noise ratio and spatial 
resolution. A sufficient spatial coherence is important to 
avoid cancellations of the contribution from different re- 
gions in the chromosphere evolving with different phases. 
If we consider the expected capabilities of the next gen- 
eration of solar telescopes (like the European Solar Tele- 
scope, EST, or the Advanced Technology solar Telescope, 
ATST), we can aim at observing the emergent Stokes 
profiles of Fig. [Jj with a 10 s cadence (upper panelj^- 
However, with the present telescopes and instrumenta- 
tion, we are forced to integrate in time and/or space to 
detect the scattering polarization signals. If we degrade 
the temporal resolution of our results to an integration 

6 Using EST (telescope diameter of 4 m, instrumental efficiency 
around 10%) and considering a spectral resolution of 30 mA, a 
spatial resolution of 0.1 arcsecs and an integration time of 1 s (ten 
times better than needed), it would be possible to observe the linear 
polarization of the 8542 A line (line to continuum ratio of ~ 0.2) 
at the level of Q/I ~ 10~ 3 with a confidence of 3a over the noise. 



time of 1 min (middle panel of Fig. [JJ and 3 min (lower 
panel of Fig. [7]) we clearly see that the time evolution 
becomes more difficult to detect. In the last case, the pro- 
files are already so smoothed that the original features 
are completely lost, both in the spectral and temporal do- 
mains. The amplitude of the integrated signals are lower 
than in the original 10 s sequence by a factor of 2 (see the 
color scales). However, integration during time intervals 
~1 min could reveal the amplification/modulation effect 
if we capture spectro-polarimetric signals similar to the 
ones showed in the middle row panel. 

4.7. The effect of a decreasing velocity on the averaged 

profiles. 

We also calculated what happens to the emergent aver- 
aged profiles in the strongly dynamic case (with 15 min- 
utes of integration, emulating an observation) when we 
gradually reduce the velocity field by a constant scaling 
factor F, keeping the rest of atmospheric magnitudes un- 
perturbed (see Fig. [8]). As expected, we find that the po- 
larization amplitudes decrease in proportion to F, from 
the original case, with F = 1, towards the static case, 
with F = 0. Note that the core of the 8498 A line goes 
through zero for a certain F value (near F = 0.6). Thus, 
depending on the magnitude of the velocity gradients, its 
linear polarization amplitude will be positive or negative. 
This fact suggests an additional way to diagnose velocity 
gradients along the line-of-sight. However, it is impor- 
tant to keep in mind that this sensitivity also depends 
on the variations in density and temperature, as shown 
in Sec. 

Furthermore, the variation of the Q/I amplitudes is 
not linear with F. The change is small for small F, is 
larger for intermediate values of F, and again becomes 
smaller for the largest F, tending to saturation. 

5. CONSIDERATIONS ON THE HANLE EFFECT 

For magnetic field diagnostics with the Hanle effect 
it is often nece ssary to know the zero-field polarizatio n 
reference (e.g., lStenflol[l99llTruiillo Bueno et al.ll2004f l. 
That is what we have tried to do in previous sections, cal- 
culating and explaining the temporal evolution of the lin- 
ear polarization profiles in chromospheric dynamic sim- 
ulations. Ideally, this reference has to be computed un- 
der the same thermodynamical and dynamical conditions 
than in the real Sun but without magnetic field. As the 
Hanle effect often depolarizes the linear polarization sig- 
nals, the difference between the observation and the zero- 
field calculation can be associated to a magnetic field by 
adjusting the magnetic field topology and intensity. The 
key point is that the reference amplitude must be as pre- 
cise as possible. If it is imprecise, variations in the Stokes 
profiles can be associated to a magnetic field when they 
were really due to uncertainties in another magnitudes, 
like the temperature or the velocity field. Due to this 
reason, the fact that the solar chromosphere is a highly 
dynamic medium brings some complications for the use 
of the Hanle effect as a diagnostic mechanism. 

A strategy to avoid the above-mentioned problem 
is known as the line ratio technique. It consists of 
finding a pair of spectral lines whose thermodynami- 
cal behaviour is identical but whose sensitivity to the 
magnetic field is different in some range of magnetic 
field intensity or inclinations (e.g., iStenflo et al.l 119981 ; 
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Fig. 9. — Black solid lines correspond to the strongly dynamic case. Red lines correspond to weakly dynamic case. Black dashed lines cor- 
respond to the strongly dynamic case but doing the calculations without micro-turbulent velocity. The line ratio £3 = (Q/I)pp 98 /(Q/I)pp 42 
is not shown because it is very similar to qi and because can be obtained from the other two ratios. 



iManso Sainz et all 12004) . In that case, the ratio be- 
tween the polarization amplitudes should only change 
due to variations in the magnetic field, thus allowing 
us to measure it after a suitable c alibra tion. As shown 
byH 

anso Sainz fc Truiillo Buenol (|2010ft . the mam mag- 
netic sensitivity difference among the lines of the Ca II 
IR triplet is between the A8498 line (which is sensitive to 
field strengths between 0.001 G and 10 G) and any of the 
A8662 and A8542 lines (which react mainly to sub-gauss 
magnetic fields and up to 10 G in the latter spectral 
line). Unfortunately, while the line-cores of the A8662 
and A8542 lines originate in similar atmospheric layers, 
the A8498 line-core originates at significantly deeper at- 
mospheric layers (see Figs. 1 and 5). Nevertheless, we 
have found useful to plot in Fig. [9] the time evolution of 
the following polarization line ratios: 
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where the super-index indicates the central wavelength of 
the transition in A. These quantities were calculated for 
each simulation considered before (weakly and strongly 
dynamic cases). The more stable they are, the more 
useful they will be for inferring the magnetic field. 

We obtain that , in average, §\ = 0.15 ± 0.10 and 
Qi = 0.16 ± 0.14 for the strongly and weakly dynamic 
cases, respectively (lower panel in Fig. [9]). The sudden 
shape variations (including maximum amplitudes pass- 
ing by zero) of the 8498 A line induce large instantaneous 
excursions on g\ . As expected, a more stable line ratio 
is obtained for the second pair of transitions, which are 
precisely the ones that originate at similar chromospheric 
heights. We find g 2 = 1.06 ± 0.11 and g 2 = 1.00 ± 0.09 
for the strongly and weakly dynamic cases, respectively 
(upper panel in Fig. If we repeat the calcula- 

tions setting to zero the micro-turbulent velocity in the 
strongly dynamic case we obtain q\ = 0.22 ± 0.09 and 



02 = 0.97 ± 0.12 (dashed black lines in Fig. EJ). These re- 
sults indicate that the g~ 2 line ratio shows a relatively sta- 
ble behavior against variations of the velocity and tem- 
perature in the solar atmosphere. Consequently, in prin- 
ciple, q 2 could be used as a suitable line ratio to estimate 
the magnetic field from spcctropolarimctric observations 
of the A8662 and A8542 lines. 

Regarding the sensitivity of these lines to the mag- 
netic field and their applicability for the diagnostic of 
magnetic fields through the Hanlc effect, several consid- 
erations have to be taken into account. First, the micro- 
turbulent velocity has a small influence on the averaged 
amplitudes and line ratios. Second, once the magnetic 
field is included in the calculations, the Hanle effect typ- 
ically operates at the line center for static cases. How- 
ever, in a dynamic situation there is not a preferred line 
center wavelength. As the maximum of the absorption 
and dispersion profiles occurs at different Dopplcr shifted 
wavelengths, the Hanlc effect will operate in a small 
bandwidth around the l ine core. Third, according to 
the st atic calculations bv IManso Sainz fc Truiillo Buenol 
(|2010f K for chromospheric magnetic fields stronger than 
0.1 G in the "quiet" Sun, the Q/I signal of the 8662 
A line is expected to be Hanle saturated. Thus, varia- 
tions between 0.1 and 10 G could be measured with g 2 , 
being produced by changes in the linear polarization of 
the 8542 A line. Unfortunately, the fluctuations we see in 
Fig. [9] (exclusively due to the dynamics) have amplitudes 
of the same order of magnitude than those expected from 
the investigations of the Hanle effect in static model at- 
mospheres (exclusively due to the magnetic field). More 
realistic results will be obtained when carrying out cal- 
culations of the Hanle effect of the Ca n IR triplet in 
dynamical model atmospheres. In any case, it is clear 
that for exploiting the polarization of these lines, we need 
instruments of high polarimetric sensitivity. 

6. CONCLUSIONS 

The results presented in this paper indicate that the 
vertical velocity gradients caused by the shock waves that 
take place at chromospheric heights in the HD models of 
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Caiisson & Stein (1997; 2002) have a significant influence 
on the computed scattering polarization profiles of the 
Ca II IR triplet. They show changes in the shape of the 
Q/I profiles of the three IR lines and clear enhancements 
in their amplitudes, as well as changes in the sign of 
the Q/I signal of the A8498 line. Interestingly enough, 
such modifications with respect to the static case are 
evident, not only in the temporally resolved Q/I profiles 
(e.g., see Fig. 2), but also in the temporally-averaged 
(Q) / (I) profiles (e.g., see Fig. 4). This is true even 
with moderate macroscopic plasma velocities, simply due 
to the presence of strong vertical velocity gradients like 
the ones produced by shock waves. This may explain 
why the above-mentioned modifications of the scattering 
polarization profiles of the Ca II IR triplet are present 
not only in the strongly dynamic simulation (Carlsson & 
Stein 1997) but also in the weakly dynamic one (Carlsson 
& Stein 2002). 

Our investigation points out that the development of 
diagnostic methods based on the Hanle effect in the Ca 
II IR triplet should take into account that the dynamical 
conditions of the solar chromosphere may have a signifi- 
cant impact on the emergent scattering polarization sig- 
nals. This complication could be alleviated through the 
application of line ratio techniques. In Sec. [5] we have 
concluded that the ratio between the polarization ampli- 
tudes of the A8542 and A8662 transitions would be the 
best line-ratio choice. However, even in the absence of 
magnetic fields, the small fluctuations we see in the value 
of such a line ratio in dynamical model atmospheres could 
be confused with the presence of magnetic fields in the 
range between 0.1 and 10 G. Further work is necessary 
at this point. 

In any case, the fact that realistic macroscopic velocity 
gradients may have a significant impact on the scattering 
polarization profiles of the Ca n IR triplet is interesting 
and important for the diagnostic of the solar chromo- 
sphere Q On the one hand, it provides a new observ- 
able for probing the dynamical conditions of the solar 
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chromosphere (e.g., by confronting observed Stokes pro- 
files with those computed in dynamical models) . On the 
other hand, the exploration of the magnetism of the quiet 
solar chromosphere via the Hanle effect in the Ca n IR 
triplet (either through the forward modeling approach or 
via foreseeable Stokes inversion approaches) would have 
to be accomplished without neglecting the possible effect 
of the atmospheric velocity gradients on the atomic level 
polarization. 

Several points are still unanswered after this work. 
First, we need to investigate the sensitivity to the 
Hanle effect of the Q/I and U/I profiles of the Ca n 
IR triplet using magnetized and dynamical atmospheric 
models. Second, we have to investigate whether our one- 
dimensional radiative transfer results remain valid af- 
ter considering realistic three-dimensional models, such 
as those r esulting from magneto-hydrodynamical simula - 
tions (e.g. lWedemever et alJ200J : Iheenaarts "et~al|[2009] ). 
We would tentatively expect that the strong stratifica- 
tion that gravity imposes in the solar atmosphere facili- 
tate shocks that propagate mainly in the vertical direc- 
tion, but with a reduced strength, given the increased 
degrees of freedom. 

Finally, we mention that our results could be of poten- 
tial interest in other astrophysical contexts. For instance, 
the mechanism of polarization enhancement due to the 
presence of shocks might well be the explanation for the 
changing amplitudes of the linear polariz ation signals 
repor ted in variable pulsating Mira stars (jFabas et al.1 
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7 The physical mechanism is general, but its observable effects lines, such as the Ca II K line and Lyo. 

are expected to be significantly less important in broader spectral 



